One-body energy dissipation in fusion reaction from mean-field theory 
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Information on dissipation in the entrance channel of heavy-ion collisions is extracted by macro- 
scopic reduction procedure of Time-Dependent Hartree-Fock theory. The method gives access to 
a fully microscopic description of the friction coefficient associated with transfer of energy from 
the relative motion towards intrinsic degrees of freedom. The reduced friction coefficient exhibits 
a universal behavior, i.e. almost independent of systems investigated, whose order of magnitude is 
comparable with the calculations based on linear response theory. Similarly to nucleus-nucleus po- 
tential, especially close to the Coulomb barrier, there are sizable dynamical effects on the magnitude 
and form factor of friction coefficient. 

PACS numbers: 25.70.Jj, 21.60.Jz 



I. INTRODUCTION 

The discovery of deep inelastic collisions in the 1970s, 
in which a large amount of kinetic energy and angular 
momentum is dissipated from the relative motion to in- 
trinsic excitations of colliding nuclei, brought us the con- 
cept of friction in nuclear physics After this discov- 
ery, Gross and Kalinowski developed the surface friction 
model (SFM) Q that takes into account this dissipative 
effect by introducing friction forces in a Newtonian equa- 
tion. More generally, dissipative aspects in self-bound 
system such as nuclei have important role on many phys- 
ical phenomena such as fission, fusion reactions and giant 
resonances. 

In most actual models that include dissipation and that 
are used to describe nuclear reactions at energies around 
the Coulomb barrier, dissipation mechanism is assumed 
to be a one-body type, where energy dissipation is caused 
by collision of nucleons with mean-field wall and by nu- 
cleon exchange between the two partners of the reac- 
tion. These effects are known as the wall-and-window 
formula [3, B| and are essentially based on classical con- 
sideration. The order of magnitude of parameters related 
to dissipation are generally adjusted to reproduce exper- 
iments. However, large uncertainty on those parameters 
exists Q. Therefore, description of energy dissipation 
from quantum microscopic models is highly desirable. 

It is know that time-dependent Hartree-Fock (TDHF) 
model [1, 0] includes one-body dissipation mechanism 
from microscopic point of view, because of the treatment 
of self-consistent mean-field. It is worth mentioning that 
so-called fusion window problem due to underestimation 
of energy dissipation in old TDHF calculations has been 
solved by including spin-orbit interactions and time-odd 
terms in the energy density functional [1, [l(| as well as 
by breaking symmetries. Now three-dimensional TDHF 
calculations including full Skyrme effective interaction, 
that are used in recent static Hartree-Fock calculations, 
are expected to provide a better description of dissipative 
aspects 

da eh mm. u P t now, there are only a few 

studies dedicated to extracting friction coefficient associ- 



ated with one-body energy dissipation from microscopic 
mean- field approach [|| [lj] . The aim of this article is to 
investigate energy dissipation mechanism in detail from 
microscopic point of view by using state of art TDHF 
theory. 

In Ref. [l5|, [l8| , a method of extracting nucleus- nucleus 
potential and friction coefficient from TDHF has been 
proposed. This method, called Dissipative Dynamics 
TDHF (DD-TDHF), assumes that mean- field evolution 
can be properly reduced to one-dimensional dissipative 
dynamics for relative distance R between nuclei. In Ref. 
[l8| |, nucleus- nucleus potential for symmetric and asym- 
metric reactions has been systematically extracted from 
the DD-TDHF. In this work, we focus on friction coeffi- 
cient, which also comes as an output of the macroscopic 
reduction procedure. 

The paper is organized as follows. In section II, we 
briefly explain the DD-TDHF method. In section III, we 
illustrate main properties of extracted friction coefficient 
for different systems. In section IV, we discuss a method 
to estimate intrinsic excitation of colliding nuclei, and 
compare the result to the dissipated energy. A summary 
is given in section V. 



II. DESCRIPTION OF THE 
DISSIPATIVE-DYNAMICS TDHF METHOD 

Dissipative-Dynamics TDHF (DD-TDHF) method, 
originally proposed in Ref. @, [l5[ , rely on the hypothesis 
that complex microscopic mean-field evolution of head- 
on collisions can be accurately reduced into a simple one- 
dimensional macroscopic evolution given by, 
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R and P denote here the relative distance and rela- 
tive momentum between two nuclei, respectively, and are 
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computed from mean-field theory according to the pro- 
cedure described in Ref. 0i 19- In the second equation, 
V(R) and j(R) denote the nucleus- nucleus potential and 
friction coefficient, respectively, while term — c£r(P 2 /2/z) 
arises from possible relative distance dependence of the 
reduced mass fi(R). The nucleus- nucleus potential for 
symmetric and asymmetric reactions deduced from DD- 
TDHF method has been systematically investigated in 
Ref. [H[ . Technical details of the method are extensively 
discussed in this reference and we recall here the main 
conclusions of the potential study: (i) We have shown 
in Ref. that the second term in the right hand side 
of Eq. ([2]), i.e., term involving i?-derivative of reduced 
mass, has a minor role on the extracted nucleus- nucleus 
potential and can be neglected. This is also the true 
here for relative distance larger than the barrier posi- 
tion, denoted by Rb- In the following, we mainly focus 
on R > Rb, where the effect of the reduced mass R- 
dependence on friction coefficient is negligible, (ii) The 
DD-TDHF method gives similar results as the density- 
constrained TDHF (DC-TDHF) method at similar cm. 
energy [j"rjl |. (iii) As beam energy increase, the extracted 
potential tends towards the frozen density (FD) approx- 
imation, which is also expected due to the fact that den- 
sity has no time to reorganize as the nuclei collide, (iv) as 
cm. energy approaches the Coulomb barrier, dynamical 
effects induce a reduction of the barrier height, which is 
in agreement with fusion threshold deduced from TDHF 
in Ref. [TtJ and is in close agreement with the experi- 
mental barrier height, see Fig. [TJ In this figure, the solid 
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FIG. 1: Barrier height extracted from DD-TDHF minus ex- 
perimental barrier height as a function of extracted barrier 
height for reactions indicated in figure. Vb is deduced from 
high cm. energy (solid line) and from low cm. energy TDHF 
(dashed line) reactions, respectively (Values of different bar- 
rier height can be found in table I of Ref. pl|). 

line corresponds to potential height extracted using high 
cm. energies TDHF trajectories (E c . m , » Vb), whereas 
the dashed line stands for barriers obtained when cm. 
energy used in DD-TDHF approaches the Coulomb bar- 
rier (-Ec.m. ~ Vb). Dynamical reduction of barrier height 
is clearly seen for all reactions. Due to this reduction, 



estimated barrier height at low energy becomes much 
closer to barriers extracted from experiments. In sum- 
mary, investigations of Ref. [l8[ have shown the simple 
macroscopic reduction method of DD-TDHF can provide 
a useful tool to infer nucleus- nucleus potentials from mi- 
croscopic mean-field theory. Not only the potential, but 
employing the same method, we can extract informa- 
tion about one-body dissipation mechanism, which is the 
topic of this paper. 

For TDHF calculations, we use three-dimensional 
TDHF code developed by P. Bonche and coworkers with 
the SLy4d Skyrme effective force [lT| . The mesh sizes in 
space and time are 0.8 fm and 0.45 fm/c, respectively. 
As TDHF initial conditions, we solve static HF equa- 
tions 0, [2(| with the same effective force and the same 
mesh size as in TDHF. The initial distances between two 
nuclei are set between 16 and 22.4 fm. We assume that 
colliding nuclei follow the Rutherford trajectory before 
they reach the initial distance for TDHF calculations. 



III. FRICTION COEFFICIENTS FROM 
MEAN-FIELD DYNAMICS 



Friction coefficient at cm. energy well above 
the Coulomb barrier 



Besides nucleus-nucleus potential, one-dimensional 
macroscopic reduction gives also access to the friction 
coefficient j(R)- We first focus on friction coefficients 
extracted from TDHF calculation when cm. energy is 
well above the Coulomb barrier, i.e. for which extracted 
potentials identify with FD potentials |18|. Possible cm. 
energy dependence of dissipation will be discuss later. 

Fig. [21 presents an important results of this work. 
The reduced friction coefficients, defined as j3(R) = 
■y(R)/[i(R), are extracted from a fully microscopic the- 
ory without any free parameter adjustment or adia- 
batic/diabatic assumption. Reduced friction coefficient 
is systematically shown for the mass symmetric reactions 
(upper panel) 40 Ca+ 40 Ca, 48 Ca+ 48 Ca and mass asym- 
metric reactions (lower panel) 16 O+ 40 ' 48 Ca, 16 O+ 208 Pb, 
40 Ca+ 48 Ca, 40 Ca+ 90 Zr, as a function of relative dis- 
tance. As the colliding nuclei approach toward each 
other, reduced friction coefficient monotonically increases 
and is almost independent of colliding system. Therefore, 
a universal behavior of the friction coefficient is observed 
in collisions at cm. energies well above the Coulomb 
barrier. 



B. Energy dependence of friction coefficient 

As cm. energy approaches the Coulomb barrier, dy- 
namical effects such as deformation and neck formation 
can take place. These aspects are automatically included 
in dynamical mean-field calculation and have been shown 
to systematically reduce the barrier compared to the high 
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FIG. 2: Reduced friction coefficient f3(R) — j(R)/fi(R) as a 
function of R divided by the Coulomb barrier radius Rb for 
mass symmetric (upper panel) and mass asymmetric (lower 
panel) reactions. 



when density has no time to reorganize in the entrance 
channel. 

On the other hand, for E c .m. — 55 and 57 MeV, which 
are close to the Coulomb barrier, friction coefficient ex- 
hibits sizable energy dependence. The lower the energy is 
the larger the magnitude of friction. Radial dependence 
of friction coefficient is very different from that at high 
energies. At low energies, friction coefficient shows a peak 
near the Coulomb barrier and decreases to smaller val- 
ues inside the Coulomb barrier, which is located around 
R = R B ~ 10 fm. 
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FIG. 3: Reduced friction coefficient as a function of R/Rb 
for the 40 Ca-|- 40 Ca reaction at different cm. energies. 



energy limit. Similarly to potential energy landscape, we 
do expect dissipation to be modified as the beam energy 
decreases. Fig. [3] shows the extracted friction coefficient 
obtained using different cm. energies between 55 and 
100 MeV for the 40 Ca+ 40 Ca reaction. We note that the 
Coulomb barrier energy of this reaction extracted at cm. 
energy E c , m , — 55 MeV is about 53 MeV. Similarly to 
potential landscape, which tends to the FD case as the 
beam energy increases, we observe that the magnitude 
of friction coefficient does not change as the cm. energy 
increases between E c . m . = 90 MeV and 100 MeV. The 
associated dissipation corresponds to limited situations 



FIG. 4: Number of nucleon transferred from projectile to tar- 
get iV£ ans [.R(t)] as a function of R for the 40 Ca+ 40 Ca reaction 
at the same cm. energies as in Fig. [3] 



The enhanced dissipation around the Coulomb barrier 
energies is partly due to early neck formation accompa- 
nied by an increase of nucleon exchange. This is illus- 
trated in Fig. [H where the number of nucleon transferred 
from one nucleus to the other is shown as a function of R 
for the 40 Ca+ 40 Ca reaction at same energies as in Fig. [3] 
For such a symmetric reaction, the separation plane is 
located at x = 0. The number of nucleon initially in the 
projectile (taken by convention at initial time at x < 0) 
having passed through the separation plane, i.e., trans- 
ferred to the target, is estimated by, 

N&JR®] = J d 3 xpp(r,t)0(x), (3) 

where pp (resp. pr) are defined through 

p P/T (r,t)= \Mr,t)\ 2 , (4) 

ieP/T 

with (f>i (r, t) being the single-particle wave function ini- 
tially in the projectile (resp. target) and propagated 
through the mean-field. 9{x) is the step function equal to 
zero for negative x and 1 for positive x. We can equiva- 
lently define number of nucleons transferred from target 
to projectile, denoted by N^. ans . For symmetric reac- 
tions, N^. ans is equal to -ZV t ^ ans at all time. For R > Rb, 
enhancement of dissipation observed in Fig. [3] is strongly 
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correlated to increase of particle exchange. Such a strong 
correlation is indeed expected from the window dissipa- 
tion mechanism in which the magnitude of friction 
coefficient is proportional to the window area, and hence 
the number of nucleons exchanged through the window. 
From strong similarity between Figs. [3] and 21 one can 
conclude that the main source of dissipation at large dis- 
tance (R > Rb) is due to nucleon exchange. 

Enhancement of dissipation at low cm. energy is sys- 
tematically observed as seen in Fig. [5J where the re- 
duced friction coefficient (3(R) is plotted as a function 
of relative distance. It also appears that enhancement 
of dissipation before the barrier R > Rb is nearly in- 
dependent of systems investigated. For instance near 
barrier R/Rb = 1, the magnitude of friction at high 
energy is about f3(R) = 0.8 x 10 21 s _1 , while at near bar- 
rier energies, its magnitude is about four times larger 
(3(R) = 3 x 10 21 s -1 , for all systems investigated. 

For all the systems investigated, calculations show that 
the magnitude of friction coefficient rapidly decreases for 
decreasing relative distances R < Rb- We believe that 
this unphysical behavior is due to simple macroscopic re- 
duction procedure, which breaks down at low energies, 
when colliding nuclei begins to overlap strongly. In this 
case, complex dynamical effects such as onset of nuclear 
deformation 18| and non-Markovian effects become im- 
portant in dissipation mechanism, which are not incorpo- 
rated in the simple reduction procedure presented here. 
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FIG. 5: Reduced friction coefficient (3(R) = j(R)/fi(R) as a 
function of R divided by Coulomb barrier radius Rb for mass 
symmetric (upper panel) and mass asymmetric (lower panel) 
reactions. In each case, cm. energy used to extract friction 
coefficient corresponds to the Coulomb barrier energy. 



C. Comparison with other models 



There is large uncertainty between microscopic and 
phenomenological description of nuclear dissipation [f|. 
Therefore it is of great interest to provide an accurate 
description of nuclear dissipation. The macroscopic re- 
duction presented in this work serves a useful insight 
towards that goal. In this section, we compare our re- 
sults with macroscopic surface friction model (SFM) 
and microscopic calculations of Adamian et al. [2l| . The 
SFM was introduced in a classical description of relative 
motion in deep- inelastic heavy-ion collisions in Ref. 0]. 
The radial friction force was parameterized as j(R) = 
K r (dVN/dR) a , where Vn is the nuclear part of the 
nucleus-nucleus potential and K r and a are parameters. 
These parameters were fitted as K r = 4x 10~ 23 s/MeV 
and a — 2. Microscopic calculations of Ref. [2l[ is 
based on the linear response theory of nuclear dissipa- 
tion [13]. These calculations take into account time evo- 
lution of single-particle occupation factors during colli- 
sion through a consistent treatment of the collective and 
intrinsic degrees of freedom. 




FIG. 6: Comparison between reduced friction coefficient de- 
duced with DD-TDHF method for the 40 Ca+ 40 Ca reaction 
at Uc.m. = 90 MeV (solid line) and the one computed by 
Adamian et al. Q for the 64 Zn+ 196 Pt reaction at E iab = 
440 MeV (filled square). The result of SFM @] for the 
40 Ca+ 40 Ca reaction is also presented by filled triangles. 



In Fig. [51 high energy (solid line) and low energy (dot- 
dashed line) DD-TDHF results are compared both with 
the SFM case (filled triangles) and an example of linear 
response theory result obtained in Ref. [21] (filled circles). 
It is observed that at R w Rb the SFM strongly over- 
estimates the magnitude of friction coefficient compared 
to the other microscopic models, while our result agrees 
at high cm. energy with the microscopic calculations of 
Ref. [2l| (filled squares). 
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IV. DISSIPATED ENERGY, NUCLEON 
EXCHANGE AND EXCITATION ENERGY IN 
NUCLEI 

In this section, we discuss the link between internal 
excitation energy of colliding nuclei and dissipated energy 
from the macroscopic degrees of freedom. In addition, we 
give further evidence for the fact that particle transfer is 
the main source of dissipation before two nuclei reach the 
Coulomb barrier. 



The dissipated energy provides a measure for trans- 
fer of energy from relative motion to internal degrees of 
freedom during the collision. Denoting by E* the inter- 
nal excitation energy and assuming that all dissipated 
energy is converted into internal excitation energy, we do 
expect E* = -Ediss- In the following, we apply a method 
to estimate directly the internal excitation energy of pro- 
jectile and target during the early stages of collision and 
show that a good agreement is obtained with the amount 
of dissipated energy. 



A. Total Energy dissipation 



B. Estimate of internal excitation energy 



From energy conservation, we can give a simple esti- 
mate of dissipated energy i?diss in the entrance channel 
as, 



p2 

— -V DD 

2/x 



(5) 



where the second term is nothing but the relative kinetic 
energy while V DD denotes the potential extracted from 
the DD-TDHF method. Fig. [7] illustrates magnitude of 
different quantities as a function of relative distance for 
the 40 Ca+ 40 Ca reaction at E c . m . = 100 MeV. According 
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FIG. 7: Evolution of potential V DD extracted from the DD- 
TDHF as a function of R (solid line) for the 40 Ca+ 40 Ca reac- 
tion at E c . m . = 100 MeV. The difference between total cm. 
energy and relative kinetic energy is also shown by dashed 
line. Difference between two curves is nothing but dissipated 
energy in the entrance channel. 



to Eq. ((2]), we can calculate dissipated energy from the 
friction coefficient 7 using the Rayleigh formula: 

E diss (R(t)) = I dt' 7 (fl(t')) [R(t')} 2 , (6) 
Jo 

where R denotes the relative velocity deduced from 
mean-field evolution. We have checked that the above 
equation gives identical result as Eq. 



1. Single nucleus case 

For a single isolated nucleus, in the mean-field ap- 
proach, the excitation energy can be estimated using 



E* 



'MF 



MFi 



(7) 



where Sj^f 



and S j^j p 



denote ground state and excited 
state mean-field energy, respectively. For small excita- 
tions, we can approximately calculate the excitation en- 
ergy according to, 



25*~5>?-Ej0x(n,-n?). 



(8) 



In this expression, e- 1 denotes single-particle energies, e F , 
with r = n,p stands for the neutron or proton Fermi 
energy, while rii and n° are the occupation factors in the 
excited and the ground state, respectively. 



2. Dinucleus case 

In the case of two colliding nuclei, we are mainly in- 
terested in the entrance channel (R > Rb) where the 
two nuclei slightly overlaps. In this di- nuclear configu- 
ration, we define the collective variables by drawing the 
separation plane as illustrated in Fig. [8] In the entrance 
channel, colliding nuclei arc weakly excited. In order 
to calculate the excitation energy of each partner using 
the lowest order perturbation expression (|8]), we need the 
occupation factors of single-particle states. We can de- 
termine the occupation factors by constructing the over- 
lap matrix of time-dependent single-particle states in the 
right side of the separation plane [23j | , 



{i\j)R= / d 3 r</>*M)^M)0(x-xo), 



(9) 



where the separation plane is at x = xq. In this expres- 
sion, 4>i(r,t) denote single particle states originating ei- 
ther from the target or projectile. The overlap matrix in 
the left side (i\j)L is defined similarly. Occupation factors 
n a associated to the left and right sides on the separation 
plane are obtained by diagonalizing the corresponding 
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FIG. 8: Schematic illustration of the Ca+ Ca reaction in 
the entrance channel at relative distance R = 9.8 fm. Differ- 
ent curves correspond to iso-contours of total nuclear density 
while vertical line denotes the separation plane. 
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FIG. 9: Occupation factors rii = {i\i)L for neutron single- 
particle states as a function of R for the 40 Ca+ 40 Ca reaction 
at E c . m . = 100 MeV. 



overlap matrix. At large relative distance, the occupation 
factors can be grouped into two classes. We first consider 
the left sub-system which initially contains the projectile. 
States with eigenvalues of close to one will corre- 

spond to single-particle states originating from projectile 
side (left side), while those with eigenvalues close to 
correspond to single-particle states originating from the 
target (right side) and which are penetrated to the left. 
Since in the entrance channel, changes of occupation fac- 
tors are small, instead of carrying out diagonalization, 
we can use first order perturbation theory to determine 
occupation factors. In first order perturbation theory, 
eigenvalues of the overlap matrices are given by 

n a « {i\i)L = nil where \i) denotes a complete basis of 
the projectile including initially unoccupied states. As 
an example, Fig. [5] shows evolution of occupation factors 
rii for six neutron single-particle states initially corre- 
sponding to lsi/2, P3/2,±3/2, Pi/2, 4/2, ±5/2, 2si/ 2 , and 
e?3/2, ±3/2 states (with the notation £jj z ) as a function of 
R for the 40 Ca+ 40 Ca reaction at E c ' m . = 100 MeV. Oc- 
cupation factors decrease monotonically as R decreases 
but remains still close to one around the Coulomb barrier 
R B « 9.8 fm. 

In order to have a simple estimate of the excita- 
tion energy, we further assume that occupation factors 
rii = (i\i)L have the shape as the Fermi-Dirac distribu- 
tion. Then, only tail of the distribution contributes to the 
excitation energy, consequently, the excitation energy of 
the left sub-system can approximately be calculated ac- 
cording to 

A P 

E* L (t) « 25>O-40(<i|i) i -n l °), (10) 

where summation runs over the states that are initially 
occupied in the projectile. Similar expression can be 
found for the excitation of the right sub-system. In 




9.5 10 10.5 11 11.5 12 
R [fm] 



FIG. 10: Comparison of internal excitation energy E* (solid 
line) with dissipated energy -Ediss from Eq. §§§ (dashed line) 
as a function of R for the 40 Ca+ 40 Ca reaction at E c . m . = 
100 MeV. Vertical arrow indicates the position of the Coulomb 
barrier radius. 



Fig. [TU1 total excitation energy E* = E* R + E* L (solid 
line) is compared with dissipated energy E^& s given by 
Eq. © (dashed line) for the 40 Ca+ 40 Ca reaction at 
-E-cm. = 100 MeV. The agreement between two quan- 
tities is quite good all the way up to the barrier distance 
R > Rb- Near equality of the excitation energy and the 
dissipated energy provides a consistency check the fric- 
tion coefficient extracted from the TDHF simulations. 



V. SUMMARY 

Using a macroscopic reduction procedure of the mean- 
field theory proposed in Ref . [1, UH, [l|[ , we extract the 
friction coefficient associated with one-body energy dis- 
sipation in the entrance channel of heavy-ion fusion re- 
actions. The magnitude and form factor of reduced fric- 
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tion coefficient have a universal property for various reac- 
tions which are investigated. Nucleus-nucleus potentials 
obtained with the same method exhibits energy depen- 
dence. In a similar manner, magnitude and form factor 
of the extracted friction coefficient depends on the beam 
energy as well. It is observed that the rate of dissipa- 
tion increases as beam energy approaches the Coulomb 
barrier. The enhancement of dissipation rate at low en- 
ergy is a consequence of early neck formation and increas- 
ing rate of particle exchange between projectile-like and 
target-like nuclei. The order of magnitude of dissipation 
deduced from TDHF is in agreement with microscopic 
calculations of friction coefficient based on the linear re- 
sponse theory. We estimate the excitation energy in the 
entrance channel by other method. Very close agreement 
found between the calculated excitation energy and the 



dissipated energy determined from the friction coefficient 
provides a further support for the validity of the extract- 
ing procedure for dissipation that we employed. 
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